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Phase Diagram of QCD at Finite Temperatures with Wilson Fermions 
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'^Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 

Phase diagram of QCD with Wilson fermions for various numbers of flavors Np is discussed. Our simulations 
mainly performed on a lattice with the temporal size iVj = 4 indicate the following: The chiral phase transition 
is of first order when 3 < Nf < 6. while it is contimions when 7V_f = 2. For the realistic case of massless n and 
d quarks and the strange quark with = 150 MeV, the phase transition is first order. The sharp transition 
in the intermediate mass region for Np = 2 at JVj = 4 observed by the MILC group disappears when an RG 
improvement is made for the pure gauge action. 



1. Introduction 

In this article I discuss the finite temperature 
phase diagram of QCD with Wilson fermions. In 
this formulation chiral symmetry is explicitly bro- 
ken by the Wilson term even for vanishing bare 
mass. The lack of chiral symmetry causes much 
conceptual and technical difficulties in numerical 
simulations and in physics interpretation of data. 
Therefore we start our discussion with the prob- 
lem of how one defines the quark mass and the 
chiral limit in sect. 2. In sect. 3 we examine two 
fundamental problems: whether the chiral limit 
of the finite temperature transition exists at all, 
and if it does exist, what is the order of the chi- 
ral transition. The strategy we take to investi- 
gate these problems is to carry out simulations 
along the line of the critical hopping parameter 
("on-iT^'' simulation). With this method we have 
identified the chiral transition points and the or- 
ders of the transition for the cases of Np = 2, 3 
and 6. Then a more realistic case of two massless 
quarks and one light quark is discussed. 

Recently the MILC group reported an unex- 
pected finding of a first-order transition at the 
intermediate quark mass region for the Np = 2 
case. Another strange phenomenon is that the 
quark mass calculated at high temperatures does 
not agree with that at low temperatures when 
f3 < 5.5. These points are discussed and possible 
resolution with an RG improved action will be 
given in sect. 4. Finally conclusions will be given 
in sect. 5. 



2. Definition of quark mass 

The chiral property of the Wilson fermion 
action was carefully investigated through Ward 
identities for axial vector currents by Bochicchio 
et al.[l]. Independently we proposed[2] to define 
the quark mass by 

2mq < I P I TT >= -TO^ < I ^4 I TT > 

where P is the pseudoscalar density and A4 the 
fourth component of the local axial vector cur- 
rent. We first showed in quenched QCD at zero 
temperature that the pion mass vanishes at the 
point which is almost identical to that where the 
quark mass vanishes. Maiani and Martinelli[3] 
also showed this for a quark mass similarly de- 
fined through the axial-vector Ward identity. 

It has been shown that the value of the quark 
mass does not depend on whether the system is 
in the high or the low temperature phase through 
simulations at /3 = 5.85 in the quenched QCD [4] 
and at /3 = 5.5 for the Np = 2 case[5] 

Fig.l shows numerical results of the critical 
value of the hopping parameter for Np = 2. 
Result for (3 < 5.0 are obtained by two methods: 
either from the vanishing point of the extrapo- 
lated TO^ or from that of the extrapolated niq in 
the confining phase [6-11]. The data at /3=6.0 
and 10.0 are new, which are determined from 
TOg=0 in the deconfining phase. A characteris- 
tic feature in this figure is a sharp change of K^. 
around f3 « 5.0 signalling a transition from the 
strong coupling region to the weak coupling re- 
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Figure 1. The line of for Np = 2. 



gion. 

3. Chiral transition 

In Fig.l a schematic line is drawn for the line 
of transition Kt(P) separating the high tempera- 
ture regime from the low temperature regime for 
a fixed temporal lattice size Nf. Whether the Kt 
line crosses the line is not trivial. As was first 
noticed by Fukugita et al.[l2], the Kt line creeps 
deep into the strong coupling region. 

To solve this problem, we take the strategy of 
performing simulations on the line KdP) starting 
from a ft with the system in the high temperature 
phase and reducing f3. We call this method "on- 
i^c" simulations. The number of iteration TVinv 
needed for the quark matrix inversion provides a 
good indicator to discriminate the two phases[13] 
because there are zero modes around Kc in the 
low temperature phase, while none exits in the 
high temperature phase. Combining this with 
measurements of Polyakov loop, plaquette and 
hadron screening masses, we identify the cross- 
ing point ftct of the KT(ft) line with the Kdft) 
line. We also check that the crossing point thus 
determined is consistent with an extrapolation of 
the line KT(ft) toward the chiral limit. From the 
behavior of physical quantities toward ftct, we are 



also able to determine the order of the chiral tran- 
sition by this method. 

The values of Kc{ft) are slightly different (at 
most of the order of 0.01) depending the definition 
of the chiral limit. They also depend on Np. We 
find that the difference due to Np is of the same 
order of magnitude as that due to the difference 
of definition. We take the largest (farthest) 
in order to see the existence of the crossing point 
since this is the most stringent condition. 

3.1. Np = 2 

Several groups[12.7.15.5.9-ll] investigated the 
problem of the finite temperature transition for 
the Np = 2 case. The results for KTift) for N = 
4 and 6 together with those for Kdft) are plotted 
in fig. 2. The curves are for guiding eyes. 

We identified the crossing point (3ai by the "on- 
ife" simulation[9] to be at ft.t ~ 3.9 - 4.0 for 
TVt = 4 and at ft^t ~ 4.0 - 4.2 for N = 6. The 
existence of the crossing point is important, be- 
caiise otherwise we can not reach the chiral limit 
in the confining phase[12]. 

The change of 7Vi„„ with (i is gradual and we 
have not observed two-state signals. Furthermore 
the pion screening mass squared gradually de- 
creases towards zero as the chiral transition point 
is approached[9]. See also Fig. 1 in [14]. These 
results imply that the transition is continuous for 
Np = 2. This is in sharp contrast with the case 
of Np = 3 and 6 where clear two-state signals are 
observed. 

For Nt = 18 with the spatial size 18^ x 24, 
we find that the transition is at ft^t 4.5 - 5.0[9]. 
Although the spatial size is not large enough, this 
result suggests that the shift of ft^t with N is 
small. 

3.2. Np =3 [14] 

The phase diagram for Np = 3 with Nf = 4 
together with the Kt line for Np = 2 is shown 
in Fig. 3. We found clear two-state signals at the 
crossing point located at ftd = 3.0 by the "on- 
Kc' simulation[9] as shown in Fig. 4. 

An interesting question is to determine the crit- 
ical value of the quark mass m"'^* up to which the 
first order phase transition persists. 

We observed clear two state signals at /3 = 4.0, 
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Figure 2. The transition/crossover lines Kt for TYj = 4 and 6 together with the Ka hne for Np = 2. 



4.5 and 4.7, while for /j = 5.0 and 5.5 no such sig- 
nals have been seen. The value of the quark mass 
m,qa at the transition of ft = 4.7 equals 0.175(2). 
Using ~ 0.8 GeV obtained from hadron spec- 
trum calculations[14], we then obtain a bound on 
the critical quark mass to"'* > 140 MeV or equiv- 
alently (to^/to^)™* > 0.873(6) 

We note that these values are much larger than 
those for the staggered fermion case for which 
the critical mass is from 10 to 40 MeV and 



0.42 - 0.58. (TOgft 
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0.075[16a8]. We have used the value of a" 
500MeV a.t ft = 5.1 - 5.2 for 7Vjr=4[19], because 
the data of a for Np = 3 are not available.) 

3.3. Np = 6 and Np >7 

We have also examined the Np = 6 case and 
identified the crossing point at ftd = 0.3[9]. where 
we see clear two-state signals (see Fig. 4). The 
small value of ftd is consistent with our finding 
that the Kt line does not cross the line at 
finite ft when Np > 7[13]. 

It should be noted that the results that the chi- 



ral transitions for Np = 3 and 6 are of first order, 
while it is continuous for Np = 2, are consistent 
with the prediction based on universality[17]. 

3.4. Np =2+1 [14] 

Now let us discuss a more realistic case of mass- 
less u and d quarks and a strange quark. When 
the mass of the strange quark changes from zero 
to infinity the order of the phase transition should 
change from first to continuous at some critical 
mass. A very important issue is whether the 
first order transition persists np to the physical 
strange quark mass. Our strategy to investigate 
this problem is similar to the "on-iiTc" simulation; 
we set K = for n and d quarks and fix the dif- 
ference between the strange quark mass and the 
u and d quark masses to some value, and make 
simulations changing the value of ft. 

We have studied two cases corresponding to the 
strange quark mass of about 150 MeV and 400 
MeV in terms of the lattice spacing estimated 
from the p meson mass (See Fig. 3 in [14].). It 
should be noted that the physical strange quark 
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Figure 3. The transition/crossover line Kt for 
Np = 3 together with that for Np = 2 and the 
Kr line. 



mass determined from = 1020Mfiy turned 
out to be rrig ~ 150_Mey[14] for our definition of 
the quark mass. 

For both cases we observed clear two state sig- 
nals, at /3 = 3.5 for 150 MeV and at /3 = 3.9 
for 400 MeV as shown in Fig. 4. Thus our re- 
sults suggest that the transition for the physical 
strange quark mass is of first order. 

This result is apparently in conflict with the 
result for staggered quarks by the Columbia 
group[18] that no transition occurs at m-aCt = 
rrida = 0.025, m^a = 0.1 (to„ = m,i ~ 12 MeV, 
m., ^ 50 MeV using the value of a from [19]), 
which suggests that a first-order phase transi- 
tion does not occur for the physical u,d and s 
quarks. Note in this connection that the results 
for Np =2 + 1 with Wilson and staggered quarks 
are both consistent with the results for degener- 
ate Np = 3 cases, respectively, that m"'* > 140 
MeV for Wilson quarks and ?r),^/'* ~ 10 - 40 MeV 
for staggered quarks. ( It is satisfied that for 
TOo<m"'* the transition is of first-order, as it 
should be.) However, consistency between Wil- 
son and staggered quarks is an issue which should 
be investigated in future. 



4. Crossover behavior at /3 « 5.0 for Np = 2 

4.1. MILC results 

The MILC group extensively investigated the 
crossover region at /3 « 5.0 for the Np = 2 case 
and found unexpected phenomena[10.11]. For 
7Vi=4, the transition is smooth both for heavy 
and light quarks. However, in the range of in- 
termediate masses corresponding to /3 ~ 5.0, the 
transition is very sharp. This is completely oppo- 
site to what was supposed to be realized. 

In Fig.5, the MILC data[10] at /3 = 5.0 to- 
gether with ours[9] are plotted. The change of 
the Polyakov loop is sharp and has a cusp. 
Further the cjuark mass in the high temperature 
phase does not agree with that in the low tem- 
perature phase and shows peculiar behavior. This 
strange behavior was also noticed by us in various 
cases (see, for example. Fig. 1 in [13]). 

The strange behavior of the quark mass for /3 < 
5.5 implies that the high temperature phase in 
this region is far from the continuum limit. On 
the other hand, the chiral transition even on a 
Nt = 18 lattice with the spatial size 18^ x 24 
occurs at /3 < 5.5. Therefore in order to obtain 
quantitative predictions for physical quantities in 
the high temperature phase we have to perform 
simulations at Nf > 18. 

Furthermore, for Nt = 6, the MILC group ob- 
served clear two state signals in the range of inter- 
mediate mass (K = 0.17 and f3 = 5.22, K = 0.18 
and (3 = 5.02, and K = 0.19 and (3 = 4.8). They 
suggested the possibility that the transition at 

= 0.19 is a bulk transition based on their nu- 
merical results. 

Now let us carefully look at Fig. 2 again. The 
Kt line for N-t=6 starts from the line at 
(3 « 4.0 - 4.2, and initially deviates from it. We 
observe, however, that the Kt line approaches 
the Kc line at /3 « 4.8 - 5.2 and K « 0.17 - 
0.19. It is very interesting to note that this is 
exactly the region where two state signals are ob- 
served. I believe that the appearance of first order 
phase transitions in the intermediate mass region 
is closely related with this unusual relation be- 
tween the Kt line and the K^ line. It is plausi- 
ble that the first order phase transition and the 
strange behavior of the quark mass illustrated in 
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Fig. 5 are lattice artifacts. 

4.2. Possible resolutions with an RG im- 
proved action [20] 

We have suggested that the problems encoun- 
tered by the MILC group at intermediate quark 
masses are lattice artifacts. There are two pos- 
sible approaches to avoid this situation: 1) Per- 
form simulations on a lattice with large N-t. How- 
ever, we have to make simulations on a lattice 
with Nf > 18, which is very difficult. 2) Use an 
RG improved action for simulations, for which 
the approach to the continuum limit should be 
smoother. 

In connection with the second possibility, we 
proposed about ten years ago an RG improved 
action for pure gauge action of the form[21] 

S',au,e = l//(co5]W(l X 1) + caI]TV(l X 2)) 

with ci = -0.331, Co = 1 - 8ci, where W(l x 1) 
and W(l X 2) are the Wilson loops of sizes 1x1 
and 1x2, and the sum is over all oriented loops. 

With this improved action we calculated the 
string tension and hadron spectrum, and inves- 
tigated the stability of instantons and the U(l) 
problem in quenched QCD [2,22]. 

The results we interpret as improvements 
may be summarized as follows. 1) The ratio 
h-jjgl f^iM = 0.488 is of order unity in contrast to 
Aj^/ Astandard = 28.81. wliicli implies that the 
bare coupling is already an improved coupling. 
2) The approach of the plaquette to the pertur- 



bative value is much faster than the case of the 
standard one-plaquette action. 3) A first-order 
phase transition found in pure SU(5) gauge the- 
ory with the standard action disappears with the 
improved action. 4) Concerning stability of in- 
stantons. the improved action is located on the 
boundary separating the stable and unstable re- 
gions in the parameter space of the action. 

With these experiences it is worthwhile to use 
the RG improved action for full QCD. There is 
also a possibility to improve the quark action. 
However we take the Wilson action for the quark 
action as a first step, because we believe that 
the effect of improvement of the gauge sector is 
much more significant than that of the quark sec- 
tor by the following reasoning. After integra- 
tion over fermion variables, the action may be 
written as a sum over various types of Wilson 
loops. The coefficients thus obtained are numeri- 
cally much smaller at simulation points compared 
with Ci = —0.331 for the improved pure gauge ac- 
tion . 

In Fig. 5 we compare the results for the stan- 
dard action at fj = 5.0 with those for the im- 
proved action at /3 = 2.0. The inverse lattice 
spacings are identical; 0.99(6) GeV for the stan- 
dard action[7] and 1.01(4) GeV for the improved 
action. In the case of the improved action, there 
is no sudden change of the quark mass in the high 
temperature region and its value agrees well with 
that in the low temperature region. The change 
of the Polyakov loop is smooth, and the pion mass 
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Figure 5. Tho data on JV^ = 4 (filled symbols) and Nf = 6 or 8 (open symbols) lattice at /3 = 5.0 with 
the standard action[9,10] and at /3 = 2.0 with the improved action. The Polyakov loops are shown in a 
range from f} = 1.8 to 2.15 for the improved action. The Kc is from rre^ = at TVt = 8. 



squared has no cusp. 

Preliminary results in the range (3 = 1.6 - 2.15 
indicate that the transition is smooth in the in- 
termediate mass region [ft ^ 2.0) and becomes 
sharper toward the chiral transition which is lo- 
cated at /3 1.5 - 1.6. (See Polyakov loops in 
Fig. 5.) Therefore we conclude that the sharp 
transition in the intermediate mass range disap- 
pears with the RG improved action at Nt = 4. 
This suggests that the first order phase transi- 
tion observed by the MILC group at Nt = 6 dis- 
appears with this action. To confirm this, we are 
planning to extend runs to a wider range of pa- 
rameters, in particular, on a lattice with N-t = 6. 



5. Conclusions 

Investigation of the phase diagram of QCD 
with Wilson fermions is conceptually and techni- 
cally difficult due to the lack of chiral symmetry. 
Therefore we first have to grasp the gross feature 
of the phase diagram. In this aspect the research 
along the boundary of the phase diagram is im- 
portant. The nontrivial boundary, the line, 
can be determined by the vanishing point of the 
quark mass defined through the axial vector Ward 
identity. Then simulations along the line as 
well as in the strong coupling limit give useful in- 
formation for the phase structure. Main results 
obtained in this way may be summarized as fol- 
lows. For the degenerate case, the chiral phase 
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transition is of first order wlicn 3 < Np < 6, 
while it is continuous when Np = 2. For the re- 
alistic case of massless u and d quarks and the 
strange quark with mq = 150 MeV, the phase 
transition is first order. 

When Np > 7, the Kt line does not cross the 
Kc line at finite /3. In the strong coupling limit, 
when the quark mass is lighter than a critical 
mass, the chiral symmetry is not spontaneously 
broken and the quark is not confined[13]. This is 
the zero temperature phase structure. I hope we 
can discuss the phase diagram at zero tempera- 
ture in detail elsewhere. 

Recent work by the MILC group has revealed 
unexpected phenomena, which occur not along 
the boundary of the phase diagram. To under- 
stand the phenomena we have to investigate the 
phase diagram in a wider parameter space of 
the action. We find that when an EG improve- 
ment is made for the pure gauge action, both the 
sharp transition in the intermediate mass region 
for Np = 2 at A''^ = 4 and peculiar behavior of 
the quark mass in the high temperature phase for 
/3 < 5.5 disappear. 

Of course, we have to confirm the conclusions 
with larger N-t and need much more work to derive 
quantitative estimates of various physical quanti- 
ties. 
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the Grant-in- Aid of Ministry of Education, Sci- 
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